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1 Introduction 

The purpose of this note is to make some remarks on the law of the solution to a heat equa- 
tion with a fractional noise in time. Let us briefly describe the context and the motivation. 
Consider a centered Gaussian field {W{t,x),t G [0,T],x G M*^) with covariance given by 

B{W{t,x)W{s,y)) = {sAt){xAy) (1) 

where for x = (xi, ..,Xd),y = {yi, ■■■■.yd) £ and where we denoted x Ay = HiLi i^i ^ Vi)- 
Also consider the stochastic partial differential equation with additive noise 

^ = \/\u + W, tG [0,T], xGM'^ (2) 
ot 2 

7X0,x = 0, XGM^ 

where the noise W is defined by ([T]). The noise W is usually refered to as a space-time 
white noise^ It is well-known (see for example the seminal paper by Dalang [3j) that the 
heat equation ([2]) admits a unique mild solution if and only if d = 1. This mild solution is 
defined as 

u{t,x)= [ [ G{t-s,x-y)W{ds,dy), t e [0,T],x eM'^ (3) 
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where the above integral is a Wiener integral with respect to the Gaussian process W (see 
e.g. [1] for details) and G is the Green kernel of the heat equation given by 

I (4.t)-/^exp(-Mf) ift>0,xEM^ 
[ if t < 0,x G M'^. 

Consequently, the process {u{t,x),t £ [0,T],x G M) is a centered Gaussian process. The 
question we address in this paper comes from the following observation: it has been proven 
in [7] for d = 1 that the covariance of the solution ([3]) satisfies the following assertion : for 
every x G M we have 

E {u{t, x)u{s, x)) = (Vt + s- V|i - si) , for every s,t G [0,r]. (5) 
v2vr ^ ^ 



This fact establishes an interesting connection between the law of the solution ([3]) and 
the so-called bifractional Brownian motion. Recall that the bifractional Brownian motion 
(i?/^''^)tg[o,r] is a centered Gaussian process, starting from zero, with covariance 

R"^^it,s) := R{t,s) = ^ {{t^^ + s'^f - \t - sr^) , s,te [0,T] (6) 

with H G (0,1) and K G (0,1]. We refer to [1] and [6] for the definition and the basic 
properties of this process. Note that, if K = 1 then B^'^ is a fractional Brownian motion 
with Hurst parameter H G (0,1). When K = 1 and H = ^ then it reduces to the stan- 
dard Brownian motion. Relation ([5]) implies that the solution {u{t,x),t G [0,T],x G M) to 
equation ([2|) is a bifractional Brownian motion with parameters H = K = ^ multiplied by 
the constant 2~^^^. This is not the only connection between the solution of the heat 



equation and fractional processes. Another interesting link has been proven in |5]. 

Our purpose is to study the law of the linear heat equation driven by a fractional noise in 

time. That is, we will consider a Gaussian field (W^{t, x), t G [0, T], x G M*^) with covariance 

E (W^(t, x)W^(., y)) = i(t2^^ + s'"" - It - s|2^)(x A y) (7) 

where we will assume throughout the paper that the Hurst parameter H is contained in the 
interval (^, 1). Let us consider the linear stochastic heat equation 

ut = Au + W^, t £ [0,T],x gR'^ (8) 

with n(.,0) = 0, where {W{t, x))^^^Q^rp^^^^^d is a centered Gaussian noise with covariance 
([7]). The solution of ([5]) can be written in mild form as 



U{t,x) 



f [ G{t-s,x-y)W''{ds,dy) (9) 
Jo 

where the above integral is a Wiener integral with respect to the noise (see e.g [T] for 

the definition) and G is given 

The following result has been proven in [1] : 
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Theorem 1 The process a^))t(:[o,T],a;eiR'' exists and satisfies 

sup E {U{t,xf) < +00 
t&[o,T],xm'^ 

if and only if d < AH. 

This implies that, in contrast to the white- noise case, we are allowed to consider the spatial 
dimension d to be 1, 2 or 3. 

We are interested in characterizing the law of the Gaussian process U{t, .) defined by ([9]) 
and to see its relation with the bifractional Brownian motion. Li other words, we try to 
understand how relation ([5]) changes when the noise becomes fractional in time. We prove 
that the solution to the heat equation with additive noise fractional in time and white in 
space can be decomposed in terms of the bifractional Brownian motion with Hurst param- 
eters H = K = ^ but in this case some additional Gaussian processes appear. Concretely, 
we show that for every x G W^, 

{U{t,x)+Y^,t e [0,T]) (^Coi?p +xf ,t G [o,r]) , 

where Co > is an explicit constant and are independent Gaussian processes, also 

independent from U, with covariance functions (jl2p and (jl3p respectively. When H = ^, an 
interesting phenomenon happens: the process vanishes while the constant Co converges 
to —-7^ when H is close to one half and the process converges to ] for H 

close to one half (see Remark [2]). In other words, in the fractional case H ^ the solution 
"keeps" half of the bifractional Brownian motion i?2'2 while the other half "spreads" into 
two parts X^ and Y^ . 

2 Analysis of the covariance of the process U{t, .) 

Suppose that s < t and denote by 

R{t,s) = B{U{t,x)U{s,x)) 
where x G M*^ is fixed. We can write, from the proof of Theorem 2.7 in [1], 

R(t,s) = aH—^ I I \u-v\'^^~'^ {{t + s) - {u + v)y2dvdu (10) 
2\/27r Jo JQ 

with aH = H{2H - 1). 

The purpose of this paragraph is to analyze the covariance (jlOp of the solution U{t, x) and 
to understand its relation with the bifractional Brownian motion. The following proposition 
gives a decomposition of the covariance function of U{t, .) in the case d ^ 2. The lines of 
the proof below will explain why the case d = 2 has to be excluded. 
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Proposition 1 Suppose d ^ 2. The covariance function R{t,s) can be decomposed as 
follows 



Rit, s) = -l=aHCaf3 (2H - 1, + 2) [{t + 



\2H- 



where Cd = 2^^; (^{^iV) Beta function defined for x,y > by /3{x,y) = Jq t^ ^{1 

t)y-^dt and 



2V27r 



anCd 



daa 



2H-2 



((t + s)_a)-|+i-((t-s) + a)-5+i 



- r dais-ayi-^' [(t-af"-^ + {t + . 
Jo 



Proof: Fix t > s. By performing the change of variables u — v = a and u + v = b with 
a + b = 2u £ {0, 2t) and 6 - a = 2w E (0, 2s) in ([JO]), we get 



R{t,s) 



1 



2V27r 
1 



aH I \a\ 

s 



\2H-2 



{2t-a)^{2s+a) 



an 



2V27T 

pt—s 

+ / ' 

^0 







aV(-a) 

2s+a 



{{t + s)- byidbda 



2H-2 



i-a) 

2s+a 



2H-2 



{{t + s) ~ by^dbda 
{{t + s)- byidbda 



+ f a"^-^ [ 

J t—s J a 



a 

2t-a 



{{t + s)- by^dbda 



By performing the change of variables a 1— )• (—a) in the first summand, we get 



R{t,s) 



1 



2V27r 

f-t-S 

+ 



an 



2H-2 



2s— a 



{{t + s) - by^^dbda 



-2H-2 







t-s 



2s+a 



{{t + s) - by^dbda 



^2H-2 



a 

2t-a 



{{t + s)- by^dbda 



Remark 1 We can notice why the case d = 2 must be treated separately in the latter 
equation. The integral with respect to db involves logarithms and it cannot lead to the 
covariance of the bifractional Brownian motion. 
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By explicitly computing the inner integrals, we obtain 



R{t,s) = 



2H-2 



2V27r 

t-s 

+ I a 



t 

2H-2 



Jo 



-i{t + s)-b)-^ 



+1 



b=2s—a 



b=a 



da 



+ 



J 

Jo 

Jt-s 



<it + s)-b)- 



+1 



1 b=2s+a 



da 



b=a 
b=2t-a 



b=a 



da 



= -^anCd \ r a2^-2((i + s) - a)-i+^da - f a^''-\{t - s) + a)-i+^da 
2v27r uo Jo 

+aHCd \ a2^-2((i + s)- ayl+'da - a^^'-^it - s) - 0)"^+^ 
Uo Jo 

f a''''-\{t + s)-arl+'da- f a^''-''{a - {t - s))-l+'da 

Jt—s Jt—s 



+OiHCd 
1 



2V27r 



t+s 



,2H-2 



((t + s)-a)~5+ido 



j't—S 



,2if-2 



((t-s)-a)-5+ida 



where 



1 



2V27r 



[/' 

Jo 



^2H-2 



{{t + s) - a)~^+^da 



2H-2 



{{t- s) + a)~2+^da 



-,2H-2i 



. rt+s 

'{a -{t- s))-i+^da - a^^-\{t + s) - a)-l+'da 
At this point, we perform the change of variable a ^ and we obtain 

r\^''-\{t + s)-ar'^+^da = {t + sf"-i C a^''-\l-ari+^da 
Jo Jo 

and in the same way, with the change of variable a 1— )■ we obtain 

r\^''-\it-s)-a)--^+^da = {t-sf'-i C a'^'-^l-ayi+^da 
Jo Jo 



(11) 



As a consequence, we obtain 



R{t, s) = auCdP ( - 1, + 2 ] [(t + s) 



2H- 



(t-s) 



2H- 



with r[^^ given by Let us further analyze the function denoted by R\ {t,s). Note 

that for every s,t £ [0, T] 

-1 



(^^) Rf\t,s) = A{t,s) + Bit,s) 



where 

A{t,s) = auCd 
and 

B{t,s) = auCd 



((t + s) - a)'2+^da 



{{t - s) +a)~t+ida 



,2H-2 



{a-{t-s))-2+^da 



t+s 



,2H-2 



{{t + s)- aY^+^da 



By the change of variables a — t = d, we can express B as 
Bit,s) = anQ ' 



J-s Jo 

anCd [ da{s - a)-^+^ [{t - af"-"^ + {t + 
Jo 



da 



\2H-21 
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and the desired conclusion is obtained. ■ 

Let us point out that the constant Cd is positive for d = 1 and negative for d = 3. This 
partially explains why we obtain different decompositions in these two cases. Thanks to the 
decomposition in Proposition [H we can prove the following proposition. 

Theorem 2 Assume d = 1 and let U be the solution of the heat equation with fractional- 
white noise Let B2'2 be a bifractional Brownian motion with parameters H = K = ^. 
Consider (^/^)t6[o,r] be a centered Gaussian process with covariance, for s,t £ [0,T] 



R^'^ {t, s) 



1 



2V27r 
1 



H- 



au I {s-af^'"^ {t + a)2 - {t - a)2 
(s - af^'^ \{t + a)~^ + (i - a)"'- 



2^/27f Jo 

and let (^■'^)t6[o,T] a centered Gaussian process with 



da 
da, 



(12) 



covariance 



Ry"(t, s) = 2^=aH [\s - a)-2 [{t + af"-^ + {t - a)^^-^] da 
2v27r Jo 

r{s - o)-5 [{t + a)2^-i - (t - a)2^^-i] da. 



2\/27f Jo 



(13) 
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Suppose that U,X^ , and are independent. Then for every x G M"^, 



{U{t,x)+Y^,t G [0,T 
where = ^anCaf^ {2H - 1, -f + 2) 



Law 



CoSf 5 +Xf ,t G [0,T 



Proof: Let us first justify that is a covariance function. Clearly, it is symmetric and 
it can be written, for every s,t £ [0, T], as 



1 



2V27r 



R^'{t,s) = Hi {t As - a)'^''-' {{tV s) + a)~^ + {{tV s) -ay^ 



sAt 



\2H-1 



H 



l[o,t](a)l[o,s](«)(i A s - af^-^ {{t + a)-^ A (s + a)"^) da 



+H l[o,t](a)l[o,,](a)(t A s - a)'^"^ (^(t - a)' 2 A (s - a)'2j da 

and both summands above are positive definite (the same argument has been used in [2], 
proof of Theorem 2.1.). Similarly, the function R^" is a covariance. If d = 1, we have 
Cd = 2 and 



R{t, s) = 2aHf3{2H - 1, -) [{t + s)5 - (t - s)^J + R^^\t, s) 



with 



1 



2V27r 



iif^(t,s) = 2aH I {s-a)'^^-'^ {t + a)^ - {t - a)^ 
Jo L 

-2aH {s - a)^ [{t + af^-^ + {t - a 
Jo 



H {s-ay^-^ {t + a)-2+{t-a)-2 



da 
2^-2] da 
da 



-2aH [ {s- a)^ [{t + 
Jo 



\2H-2 



+ (*-«: 



,2H-2-\ 



da 



where we used integration by parts in the first integral. ■ 
In the case d = 3 we have the following. 

Theorem 3 Assume d = 3. Let B2'2 be a bifractional Brownian motion with H = K = ^ 
and let (Z/^)(g[o^r] be a centered Gaussian process with covariance R^^\t,s). Then 

{u{t,x) + CoB-2^-2,te [o,r]) {zf',te [o,r]), 

with Co defined in TheoremlM 
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Proof: We have C-? = —2. In this case we can write 



with 



R{t, s) + 2aHp{2H - 1, -) ^{t + s 



,2H- 



it-s) 



2H- 



-1 



Rf\t,s) 



Jo 



2aH I {s-af"'"^ {t + ay^ -{t-a) 



-2aH I (s — a] 2 
'o 



da 



\2H-2 



+ {t-a 



,2H-2^ 



(3) 

Note that i?} is a covariance function because it is the sum of two covariance functions. 



Remark 2 Let us understand what happens with the decompositions in Theorems andl^ 
when H is close to |. We focus on the case d = 1. The phenomenon is interesting. We 
first notice that the process vanishes in this case. The covariance of the process 
becomes 



R''^{t,s) 



1 



4\/27f Jo 



it 



it 



da 



1 



2V27r 



The constant Co = -^=aH(3{2H — 1, |) is not defined for H = ^ because of the pres- 
ence of 2H — 1 in the argument of the beta function. But the following happens: since 

is an approximation of the unity, it follows that aHf^{2H — 1, |) 



\2H-2 



2a/fl(o,i)('"')(l - u 
converges to ^ when H goes to ^. Therefore Co becomes 



1 



2V2iT 



. Therefore we retrieve the 



result in and recalled in relation In other words, in the fractional case H ^\the so- 
lution "keeps" half of the bifractional Brownian motion B2'2 while the other half "spreads" 
into two parts. 



3 Analysis of the processes and 

In this paragraph, we will prove some properties of the processes and appearing 
in Theorem [21 It would be nice to find that these processes are more regular than the 
bifractional Brownian motion B2'2 and therefore reduce the study of the solution U to the 
study oi B2'2 whose properties are relatively well-known. 
We start with the following remark. 



Proposition 2 The processes (X/^)t>o, (^t^)t>o o'^t? (Z/^)(>o from Theorems IE o-nd\^ 
2H — i self-similar. 



are 
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Proof: Indeed, for every t > s and for every c > 

\2H-1 



1 \ -1 



R^" {ct, cs) = H I {cs - af^~^ {ct + a)~^ + {ct - a)'^ 
2\/27ry Jo 



da 



2V27r 



R^" {t, s) 



where we used the change of variables a' = |. The same calculation can be done for 
and Z". U 

Let us now estimate the norm of the increments of the process in order to study 
the regularity of its paths. 

Proposition 3 For t > s we have 

E - |2) < C\t - s\-2 and E (1^^^ - Y.^p) < C\t - 

with C a strictly positive constant. Consequently, and each admit a version that is 
Holder continuous of order 5 G (0, 

Proof: Using the expression of the covariance of the process X^ , 
E|Xf-Xf|2 = R''"{t,t)-2R''"{t,s) + R''"{s,s) 



H / (t-a)'^-i {t + ay2 +{t-ay 
Jo 

+H [ {s- af"-^ \{s + a)-^ + (s - a)-^ 
Jo L 



da 



-2H {s-af^^^ (t + a)-2 +(i-a) 

Jo L 
T^{t,s) + T2{t,s)+T^{t,s) 



da 

da 



where 



{t - af^~^ [(i + a)-^ + {t- a)-^ 



da, 



and 



T2it,s) = H / {(t 
Jo 



\2H-1 



- (s - af^-^) \{t + a)-^ + {t- a)-3 



da 



Ts{t, s) = H [ {s- af"'^ \{s + a)"5 + _ a)-| _ + - [t - a)"^ 
Jo ^ 

We claim that, for every s G [0, T], t > s, we have 



da. 



Ci\t-s\^"-^ <Ti{t,s) < C2\t-sf"-'2, 
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T2{t,s) < Cslt - sf^~^ and Ts{t,s) < C4\t - s\^ 
with Ci, C2, C3, C4 strictly positive constants. Indeed, since {t + a)~2 < [t — a)~2 we obtain 

^* ^ 2H 



Ti{t,s) {t-af"~^{t-a)-^ 

For the lower bound, since {t — a)"^ > {t — s)~^, 

Ti{t,s)>H{t-s)-'2j^{t- 
Concerning the summands T2 and T3 we can write 



■2 da 



2H - 



a^^'-'da 



H 



2H-1 



{t-s 



T2{t, s)<{t- sf^~^H I {t- a)-^da = ^{t-s) 

Jo 2 



2H-1 



and 



da 



T3{t,s) < HT^"-^ [ \{s + a)~^ + {s - a)--^ - {t + a)--^ - {t - a)~^ 
Jo L 

= f ((2s)^-(t + s)5 + (i-s)i) 

The calculations are similar for . The Holder continuity of the paths follows by using 
the Kolmogorov criterium. ■ 



Remark 3 A detailed study of the term denoted by above shows that its Holder order 
cannot be improved. Therefore the regularity of the paths of and is the same as 
those of the bifractional Brownian motion B^'^ . Recall that B^'^ is Holder continuous of 
order strictly less than |. 
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